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Abstract. Solutions of the two-body problem with the simultaneous action of the solar 
electromagnetic radiation in the form of the Poynting-Robertson effect are discussed. 
Special attention is devoted to pseudo-circular orbits and terminal values of osculating 
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(— I ' elements. The obtained results complete those of Klacka and Kaufmannova (1992) and 

^ Brcitcr and Jackson (1998). 
O 

^ Terminal values of osculating elements presented in Breiter and Jackson (1998) are 
of no physical sense due to the fact that relativistic equation of motion containing only 
^ first order of v/c was used in the paper. 
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1. Introduction 

Breiter and Jackson (1998; BJ-paper in the following text) have presented analytical 
mathematical solutions of two-body problems with drag. The Poynting-Robertson effect 
(P-R effect) is one of the special cases discussed in the BJ-paper. 

The conclusion of the BJ-paper concerns the fact that analytical solution in terms of 
special functions exists for the P-R effect. However, this result is useful only in the two 
limiting cases: pseudo-circular orbits and terminal values of osculating elements. As for 
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the general case, there is a complication with calculation of infinite functional series and 
even with the divergence of the series. 

The problem of the divergence of the infinite functional scries can be overcome by 
numerical integration of the equation(s) of motion(s). This can be easily done at the 
present epoch of computers. However, it is nice, indeed, when a man can compare his 
numerical calculations with analytical solutions in the limiting cases - analytical solutions 
are very interesting. 

The aim of this paper is to discuss analytical solutions for the P-R effect in the two 
limiting cases: pseudo-circular orbits and terminal values of osculating elements. It is 
shown that the real results do not completely correspond to the results presented in the 
BJ-paper. General analytical solutions for pseudo- circular orbits presented in our paper 
are presented. 

2. Overview of the correctly presented results 

The first correct special result for the P-R effect was presented by Robertson (1937; cited 
in accordance with BJ-paper). The complete formula for the P-R effect is correctly de- 
rived in Klacka (1992a). Other, more simple correct derivations may be foimd in Klacka's 
papers: 1992b, 1993a, 1993b. These papers present also arguments as for physical incor- 
rectness in papers cited in the BJ-paper (except of Robertson). 

As for the solutions of the equation of motion for the P-R effect, we refer also to 
Klacka (1992c; error is in Eq. (10) of the paper - the right-hand side of Eq. (10) must 
contain fj, {1 — (3) instead of /x when used in the section 3), Klacka (1993c, 1993d), 
Klacka and Kaufmannova (1992, 1993 - typewriting error in Eq. (1)) (see also Klacka 
1994a). As for other papers of the author, dealing mainly with general interaction of the 
electromagnetic radiation (of the Sun; and, also, solar corpuscular radiation), we refer to: 
Klacka (1993f, 1993g (some numerical errors which may be easily found are in the last 
section; moreover, real particle should rotate around one axis - axis of rotation), 1994b), 
Klacka and Kocifaj (1994). 

3. P-R effect and analytical solutions 

At first, we put Eq. (18) in BJ-paper into a correct form: 



r = — X ivy + xy) 
at 



(1) 
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In the case of the P-R effect we have f = 1, fj, must be substituted by /xq (1 — /?) 

and at = P Ho I c: 

r = c — - — X [y + xy), 
at = /3 Ho / c . (2) 

3.1. Initial conditions 

Let the initial orbit is given by rj„, fin, hin and = 0. Eqs. (5) in BJ-paper and Eqs. 
(2) of our paper yield 

- h 
P Mo 

- ^ r- 3 

/iO (1 - /3) Tin 

J. - ^ r? 1 r- * C^) 

{ C 1 - P flQ (1 - P) Tin J 

These equations complete Eqs. (16) in BJ-paper. 



4. Special types of pseudo-circular orbits 

4.1. Increasing eccentricity 

Due to the fact that Bessel's functions Ji{x), Yi{x) can be expressed as linear combina- 
tions of cos(a; — 3 tt / 4), sin(a; — 3 tt / 4), Eqs. (16) in BJ-paper enable to fulfil the 
conditions A = B = 0. These conditions are fulfiled by conditions 

yl^ = S[{xin) . (4) 

Eqs. (3) and (4) lead to special typo of pseudo-circular orbits, as it is discussed in 
section 5 in BJ-papcr. The procedure goes in the way that giving hm we can calculate 
Xin = hin c /{(} Mo); the first of Eqs. (4) and second of Eqs. (3) yield 

^- = r^l!{^^M"^ (5) 

the second of Eqs. (4) and the third of Eqs. (3) yield 
1 _ /3 r J.2 





Transversal component of velocity may be calculated from the relation vt 

hin / f^in- 



The case e w 2/a; (Eq. (28) in BJ-paper), together with ^/xq (1 — P) P = x f3 jiq / c, 
leads to 
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(2 Pf MO 
1 - /3 c2 ■ 



(7) 



In any case the conclusion about the increasing eccentricity is not completely correct. 

The eccentricity begin to oscilatc at some state of the orbital evolution, although its mean 
value is an increasing function of time. The conclusion "The general conclusion is that 
in the pseudo-circular solution the osculating eccentricity grows systematically from the 
small but nonzero value of order at/h towards e = 1." (BJ-paper, section 5) is incorrect 
(see also section 4 of our paper). 

4-2. Oscillating eccentricity 

Pseudo-circular orbit discussed in BJ-paper is only a very special case of possible pseudo- 
circular orbits. We complete the case by all the other possibilities, where also zero value 
of eccentricity is possible. 

Let us consider a situation when mcteoroid is ejected from comet at cometary aphelion 
with zero ejection velocity. If the comet's orbital elements arc ag, eg and metcoroid's (3 
is given by the condition /3 = eq, then Eqs. (30) and (31) in Klacka (1992c) and Eq. (17) 
in Klacka (1993e; Eq. (17) must contain /xq (1 — /3) instead of fj., now) yield 



(If mcteoroid is ejected with nonzero velocity, one can use equations of Gajdosfk and 
Klacka (1999).) Then, we have 

rin = ro = ao (1 + eo) , 



If we set i?o = 0, at = /? Mo / c, equations of sections 2, 3 and 4 of BJ-paper yield 



din = ao (1 + eo) , 




(8) 



Tin = ro = ■ 



(9) 




yo = {MO (1 - /3) Tin Xo} ^ al^XQ 



(10) 



Eq. (2) yields for rj„ = 



0: y' = — y / X, and, thus 



y, 



'in = y'o = -yO / X0 = - Xq 



-4 



(11) 
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This case should correspond to the case of Klacka and Kaufmannova (1992, 1993). We 
do not treat it here, since general pseudo-circular orbit will be discussed in the following 
section. 

5. General pseudo-circular orbit 

In order of generalizing of Eqs. (10) and (11), we use Eqs. (5) of BJ-paper: 



y^Mo (1 - 13) a,„ (1 - e2J , 



P Mo 

^ _ P"^ no 1 + Bin COS /in _ 1 + e,„ COS fjn _i_ 

Vin - yo - ^2 (1 _ ^) (1 _ e2j - .3 ^ .3 i^^) 



for pseudo-circular orbits. Eq. (3) yields 

, j y/1 - P c ■ f A -4 _ k 

yo ~ \ ^ / , sin/i„ - 1 V a;.„ = - -J (13) 

for pseudo-circular orbits; u'q = — 1 / Xq for rj„ = (Eq. (11)), = — 3 / Xq + ... 
for A = S = (Eq. (4)). 

The second of Eqs. (12) and Eq. (13) yield 



+ ,/7j2 x^'/' ^1±1± eos(a;o - 3 tt / 4) 

o Xo 



Aw- Vtt / 2 Xo (3 - fc) sin(a;o - 3 tt / 4) + 

-7/2 
I 

B w / 2 a;o (3 - k) cos(a;o - 3 tt / 4) 

+ v^;r72 Xo ^ sin(xo - 3 TT / 4) . (14) 

o Xo 

If X = /i — t?, then 

Zi(x) = A Ji(x) -I- B yi(x) w xo^^^ x-i/2 (3 - k) sin(x - Xo) -|- 
-7/2 -1/2 [ 9 - 3 fc 55 + 3fc\ , , 

S'i(x) w x"^ - 8x-^ , 

Zo(x) = A Jo(x) -h S Yo{x) « Xq (3 - fc) cos(x - xo) + 

-7/2 -1/2/3 - fc 55 + 3fcl 
Xo ^ X / sin(x - xo), 

,S;(x) w - 3x-'^ + 40 x"*^ . (15) 
One obtains, then 

e cos/ = x^ (Zi(x) -h S'i(x)) - 1 w 

~ (3 — fc) Xq ''^^ a;^^^ sin(x — xo) — 8 x~^ -|- cos(x — xq) x 

.-7/2 „3/2 , 55 -h 3 fc -9/2 5/2 

(-0 J- -r „ j-o •'' 
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e sin/ = {Zq{x) + x''^ Si{x) + S[{x)) « 

« (3 — fc) Xq"^^"^ x^/^ cos(x — xo) — 2 x^"^ + sin(x — xq) x 
3 - ^-7/2 ^3/2 _ 55 + 3 A; ^-9/2 ^5/2 



{ 

4 a;-2 + (3 - fc)^ Xq ^ a;^ - 4 (3 - fc) x^"^'"^ x^'^ cos(a; - xq) 

+ x^/^ sin(x — xo) X 

r 33 (3 - fc) 55 + 3 A; _i 

< ^ — ^ Xq X + (3 — k) cos(x — XojX 

15 — 5 A; . x7/2 55 + 3 fc , _i .9/2 

- (Xq Xj - - (Xq xj 



^^^^ 

The equation for represents general type of evolution of eccentricity for pseudo- circular 
orbit in comparison with Eq. (28) in BJ-paper and papers by Klacka and Kaufmannova 
(1992, 1993). 

5.1. Time evolution of eccentricity and true anomaly 

Eqs. (16) enable us to investigate time evolution of osculating eccentricity and true 
anomaly. Wc must bear in mind that the term 'increasing time' corresponds to the term 
'decreasing x'. 

5.1.1. Extremes of eccentricity 

The local extremes of eccentricity as a function of time are given by condition (very large 
values of xq are considered, theoretically xq — > 00; A; 7^ 3) 

de^ 



dx 



= 



sin(x — Xq) « . 



(17) 



dx^ 



As for the second derivative, we have 
4 (3 - A) Xq x^/2 cos(x - xo) 



(18) 



The case A; = 3 yields no local extreme for the leading terms. Eqs. (16) yield for 
evolution of true anomaly /: 3 tt / 2 ^ tt and the value tt never occurs. 

The case A = 1 cannot be treated analytically, as it will be shown later on. 



5.1.2. X = Xo 

Putting X = Xo into Eqs. (16), (17) and (18), one obtains that the point x = xo corre- 
sponds to: i) local minimum for the case A < 3, ii) local maximum for the case A: > 3, for 
the function e{t), or, e(x). 
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If we want to find time evolution of true anomaly, we have to use Eqs. (16). The time 
t = to + At, for small positive At, corresponds to x = xo — Ax, where Ax is a small 
positive quantity (in radians). Eqs. (16) yield, then 

1 k 

e cos / « — (3 — k) — , e sin / w . (19) 

Xo Xo 

Eqs. (19) imply, on the basis of e > and |e sin/| » |e cos/|: 
lim_ / = 77 / 2 , fc < 1 , 

hm_ / = 3 TT / 2 , fc > 1 . (20) 
As for the value of eccentricity, the last of Eqs. (16) yields 

e\x = xo) « (1 - kf xf . (21) 

5.1.3. X = Xo — TT 

Putting X = Xo — -K into Eqs. (16), (17) and (18), one obtains that the point x = xq — 'k 
corresponds to: i) local maximum for the case A; < 3, ii) local minimum for the case k > 
3, for the function e{t), or, e{x). The value of the osculating eccentricity is 

e^ix = a;o - tt) « (5 - kf Xg^ + tt {5 k - 11) (5 - k) Xq^ . (22) 

If we want to find time evolution of true anomaly, we have to use Eqs. (16). 

The situation shortly before the extreme corresponds to x = xo — tt + Aa;, where 
Aa; is a small positive quantity (in radians) shortly before the extreme, and, Aa; is a small 
negative quantity (in radians) shortly after the extreme. Eqs. (16) yield, then 

e cos / Rj — (3 — k) , e sin / w — . (23) 

Xo Xo 

Eqs. (23) imply 

lim/ = 37r/2 , k<5 , 

lim/ = 7r/2 ,fc>5. (24) 

X—>Xo — TT 

Eqs. (24) yield that true anomaly is a continuous function at a; = a;o — tt - at the 
first local maximum of osculating eccentricity. 

5.1.4. a; = a;o — 2 tt 

Putting X = xq — 2 TT into Eqs. (16), (17) and (18), one obtains that the point x = 
Xo — 2 TT corresponds to: i) local minimum for the case k < 3, ii) local maximum for the 
case A: > 3, for the function e{t), or, e{x). As for the value of eccentricity, the last of Eqs. 
(16) yields 
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e^{x = a;o - 2 tt) (1 - kf x-^"^ + 2 tt (5 A; - 11) (5 - fc) x^^ . (25) 

Comparison of Eqs. (21) and (25) leads to the conclusion 

e(a:o) = e(x^ — 2 tt) , fc = ll/5, 5 
e(xo) < e(a;o - 2 tt) , 11/5 < A; < 5 

e(a;o) > e(a;o - 2 tt) , A; < 1, 1< fc < 11/5, 5 < fc (26) 

Eq. (25) shows that the case A; = 1 leads to inconsistencies: e{x — Xq — 2 tt) = - only 

the leading term can be considered. 

If we want to find time evolution of true anomaly, we have to use Eqs. (16). 

The situation near extreme corresponds to a; = xq — 2 tt + Ax, where A.t is a small 
positive quantity (in radians) shortly before the extreme, and, Ax is a small negative 
quantity (in radians) shortly after the extreme. Eqs. (16) yield, then 

e cos/ « (3 - A) — , e sin/ « ^ . (27) 

Eqs. (27) imply, on the basis of e > and |e sin/| |e cos/|: 

lim f = -K I 1 , k<l , 

x^xq — 2 tt 

lim / = 37r/2 , A>1. (28) 

X—^Xo — 2 TT 

5.1.5. X = xo — 3 TT 

Putting x = xq — 3 TT into the last of Eqs. (16) one obtains that the point x = xq — 3 tt 
corresponds to: i) local maximum for the case A; < 3, ii) local minimum for the case k > 
3, for the function e{t), or, e(x). The value of the osculating eccentricity is 

e^ix = Xo - 3 tt) « (5 - A;)2 Xq 2 + 3 TT (5 A: - 11) (5 - A;) Xq ^ . (29) 

Comparison of Eqs. (22) and (29) leads to the conclusion 

e(xo — tt) ^ e(xo — 3 tt) , k ^ 11/5, 5 
e(xo — tt) < e(xo — 3 tt) , 11/5 < A < 5 

e{xo - tt) > e(xo - 3 tt) , A; e (11/5,5)'. (30) 
As for time evolution of true anomaly, the results are analogous to Eqs. (23) - (24). 

5.1.6. Discussion 

The values of eccentricities are collected in Eqs. (26) and (30). The consequence of these 
equations is that mean values of eccentricity during the corresponding periods exhibit 
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similar properties. The cases k = 11 /5, k = 5 yield the constant values of mean eccen- 
tricities during a long time evolution. 

The case k — 1 was treated in detail in Klacka and Kaufmannova (1992). The main 
results arc presented in Figs. 1, 3 and 8 in Klacka and Kaufmannova (1992). 

We can collect the results for the cases k — 1 and fe = 5 in the statement that true 
anomaly is a discontinuous function: i) at x = xq — 2 tt m € N for A; = 1, i) at 
X = xo — -K {2 m — 1), m e iV for fc = 5. We collect the results: 

A; = 1 : lim / = TT , lim / = 3 tt / 2 , 

k = 1: lim / = 2 TT , lim / = tt , (31) 

X— >(xo — 2 7r)+ X— >(xo — 2 7r)~ 



A: = 5:lim/ = 37r/2, hm f = 2 n , 



+ 



k = 5: lim / = tt , lim / = 3 tt / 2 . (32) 

X— *(xo — 7r)~ X— >(xo — 2 tt) 

6. Terminal values of osculating elements 

As for the terminal values of osculating elements presented in BJ-paper, the results may 
be collected in two important statements: 

lim f = TT ; lim e = 1 . (33) 

x^0+ x^0+ ^ ^ 

Let us calculate other important quantities. The results are: 
lim r = ; lim a = , 

x^0+ x^0+ 

lim vt = ; lim VR = - ^ ^ ^ ^ , 

lim = ; lim £; = - oo , (34) 

x^0+ x^0+ 

where r is particle's distance from the central point mass, a - semimajor axis, vt - 
transversal component of the velocity vector, vr = r - radial component of the velocity 
vector, H - angular momentum, E - total energy of the particle with respect to the 
central point mass. 

The obtained results presented by Eqs. (33) and (34) yield important inconsistencies. 
The osculating trajectory is parabola (e = 1), the particle is situated at apoccnter (/ = tt) 
and the total energy is = — oo. Normal result is that E = for the case e = 1. So, 
there is something wrong with the results, it seems. 

The results given by Eqs. (33) and (34) are correct as for mathematical point of view 
- mathematical solutions of the discussed limits, based on the mathematical solution 
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presented in BJ-paper. However, we are not interested in mathematics as the main theme. 
We are interested in physics. Thus, the important inconsistencies presented by Eqs. (33) 
and (34) should show that physics is not completely correct. Really, physics is incorrect. 

It is a physical nonsense when a particle losses an unlimited energy within a finite time. 
The results presented by Eqs. (33) and (34) correspond to this nonphysical situation. A 
particle spirals toward r = in a finite time and its potential energy decreases in an 
unlimited value. 

Eqs. (34) yield a hint how to put the discussed inconsistencies into a correct physics. 
Since 

for < /3 < 1/3, we have to use complete form of the P-R effect - relativistic effect 
(Klacka 1992a, Eq. (140)). 

Eq. (35) yields that the form of the equation of motion containing only first order of 
v/c could be acceptable only forO<l — /J<Cl - only in this case the requirement 
t; <C c holds. However, the third Kepler's law yields = 4 tt^ {^^[3=Q (1 - /?)}~^ 
and Imifj^i- T = oo. Thus, the situation 0<1 — /3<Clis not physically interesting. 
This situation is evident also from Eq. (30) in Klacka (1992c): lim^^^- tin > 1 - no 
inspiralling toward the center occurs. 

The conclusion of this section states that the physics used in BJ-paper is not compe- 
tent to say something about the terminal values of osculating elements. Any comparison of 
the statements for various initial conditions can be done only for r final '> Tg = 2GM/c'^, 
V <^ c. (We refer also to Klacka 1994c.) As an example we may mention the time of 
inspiralling toward the point mass center, based on approximations in first order in v/c 
- the ',time of inspiralling" corresponds to r final > 200 km for the mass of the Sun. As 
for the averaged equations for osculating elements, the condition (104) in Klacka 1992d 
must be fulfilled). 

7. Conclusions 

We have completed analytical formulae (expansions) for pseudo-circular orbits for the 
Poynting- Robertson effect. 

The statement 'The zero values of our arbitrary constants do not imply that the 
osculating e = 0, but they match the case of the zero mean eccentricity (Breiter and 
Jackson 1998, section 5 on page 240) is incorect - it is not possible that any continuous 
non-negative quantity has zero mean value, unless it is identical zero. 
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Analytical results for more general pseudo-circular orbits than those discussed in 
(Breiter and Jackson 1998), were obtained in this paper. 

Wc have shown that terminal vahics of osculating elements lead to serious inconsis- 
tencies caused by the fact that relativistic equations of motion only in first order in v jc 
were used in Breiter and Jackson (1998). 

Finally, we have to stress several facts for elliptical orbits for the P-R effect containing 
only first order in vjc and in the zone of its applicability in two-body problem. 
At first, the evident result is that osculating semi-major axis is still a decreasing function 
of time - energy decreases (see Eq. (22) in Klacka 1992c). 
As for osculating eccentricity, it: 

i) may be an increasing function of time for an initial long-time interval if Eqs. (3) and 
(4) are fulfiled; 

ii) still alternates in an increasing and a decreasing functions on short- time intervals. 
As for mean eccentricity, it may be: 

i) an increasing function of time (for a suitable initial conditions: see Fig. 3 in Klacka and 
Kaufmannova (1992), or, Eqs. (3) and (4) - fc G (11 / 5,5) for pseudocircular orbits); 

ii) a constant function of time (see Fig. 4 in Klacka and Kaufmannova (1992); moreover, 
the cases fc=ll/5, fc=5 for pseudo-circular orbits); 

iii) a decreasing function of time for some initial time interval followed by an increasing 
function of time (Klacka and Kaufmannova (1992) -fee (11 / 5,5)' for pseudo-circular 
orbits); 

iv) a decreasing function of time (Wyatt and Whipple 1950 - without derivation; correct 
derivation Klacka 1992c) for a long-time interval. 

(One must be careful which type of osculating elements is used. The case ii) in mean 
eccentricity may corresponds to osculating elements defined by value /xq ~ osculating 
elements I in Klacka 1992c, or osculating elements in Klacka and Kaufmannova (1992). 
The cases i), iii), iv) and pseudo-circular orbits ii) in mean eccentricity correspond to 
osculating elements defined by value /Uq (1 ~ &) ~ osculating elements II in Klacka 
1992c, or, 'non-osculating' elements in Klacka and Kaufmannova (1992).) 
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